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(N . Abstract 

In the weak-coupling limit, kq — > oo, the partition function of simplicial 
quantum gravity is dominated by an ensemble of triangulations with the ratio 
Nq/Nd close to the upper kinematic limit. For a combinatorial triangulation 
of the D-sphere this limit is 1/D. Defining an ensemble of maximal triangula- 
' tions, i.e. triangulations that have the maximal possible number of vertices for 

a given volume, we investigate the properties of this ensemble in three dimen- 
sions using both Monte Carlo simulations and a strong-coupling expansion of 
the partition function, both for pure simplicial gravity and a with a suitable 
modified measure. For the latter we observe a continuous phase transition to a 
crinkled phase and we investigate the fractal properties of this phase. 



<D 
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■ 1 Introduction 

?— i 

Discretized models of D-dimensional Euclidean quantum gravity, known as simpli- 
cial gravity or dynamical triangulations, have been extensively studied in the past 
decade. The two dimensional model was proposed quite some time ago in Ref. Q, 
and generalized to three and four dimensions in Ref. |2|, |3| . The model is defined by 
the grand-canonical partition function 

Z(ja, K ) = ^2e ~^ Nd Z(k , N d ) (1) 

where and kq are the discrete cosmological and inverse Newton's constants, Ne> 
is the number of D-dimensional simplices, and Z(kq, No) is the canonical partition 

1 Permanent address: Institute of Comp. Science Jagellonian University, 30-072 Krakow, Poland 
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function 

Z(«o, N D ) = E e KoiY ° WbC^o, Afc) • (2) 
The micro-canonical partition function, is defined as 

W D (N ,N D )= E ^=^T E 0) 

TeT{N ,N D ) u LeC(N ,N D ) 

where the sum is over an ensemble of combinatorial triangulations T with a spherical 
topology and with iVo vertices and Njj D-simplices. For D > 2 these two numbers 
are independent. The symmetry factor Ct of a triangulation T is given by the 
number of equivalent labeling of the vertices. Including this factor is equivalent to 
summing over all labeled triangulations C(Nq, Njj). 

As a statistical system simplicial gravity displays a wealth of intriguing fea- 
tures. In D > 2 it exhibits a geometric phase transition separating a strong-coupling 
crumpled phase, where the internal geometry of the triangulations collapses, from 
a weak-coupling elongated phase characterized by triangulations consisting of "bub- 
bles" glued together into a tree-like structure or branched polymers 0]. Both in 
three |J and four dimensions || this transition is discontinuous. 

In the absence of a continuous phase transition with a divergent correlation 
length any sensible continuum limit where a theory of gravity might be defined ap- 
pears to be ruled out. This has led to the studies of modified models of simplicial 
gravity, e.g. by changing the relative weight of the triangulations with a modified 
measure ||, or by coupling matter fields to the theory [fjj. Suitable modified the 
phase structure of the model changes; the polymerization of the geometry in the 
weak-coupling limit — the elongated phase — is suppressed and a new crinkled 
phase appears ||. The crinkled phase is separated from the other two phases by 
lines of either a soft continuous phase transitions, or possible a cross-over (a diver- 
gent specific heat is not observed). Qualitatively the same phase structure, shown 
schematically in Figure 1, is observed both in three and four dimensions. 

Both in the elongated and in the crinkled phase the canonical partition func- 
tion is seen to have the asymptotic behavior 

Z( Ko ,N D ) ~ e^ N D N2~ 3 , (4) 

which defines the string susceptibility exponent 7. In the crumpled phase, on the 
other hand, the sub-leading behavior is exponential, corresponding to 7 = —00. 

Many of the properties of the model Eq. (|]]), in the weak-coupling regime 
kq > Kq, are reflected by triangulations close to the upper kinematic bound 
Nq/Nd = n™ ax . This is evident both from Monte Carlo (MC) simulations and 
from the analysis of a strong-coupling expansion (SCE) of the partition function. 
For combinatorial triangulations, where each simplex is uniquely defined by a set 
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Figure 1: A schematic phase diagram of modified simplicial quantum gravity for D > 2. 
The solid line indicates discontinuous phase transitions, the dashed lines are either soft 
continuous transitions or a cross-over. 



of (D + l) distinct vertices, the upper kinematic bound is n™ 11 * = 1/D. Although 
this bound is only saturated in the thermodynamic limit, — * oo, for each finite 
volume there exist triangulations with a well-defined maximal number of vertexes 
N™ ax (Nrj). The ensemble of these triangulations, the maximal ensemble T max , de- 
fines the micro-canonical partition function 

W max {N D ) = W D (Nr x (N D ),N D ) = cT" ( 5 ) 

In this paper we investigate the properties of the maximal ensemble, mostly 
restricted to three dimensions, using both MC simulations and a SCE, i.e. a direct 
enumeration of the micro-canonical partition function Eq. (||) . This investigation is 
separated into two parts. In the first part we study the pure (un-modified) ensemble 
and show that W max (N^) effectively separates into three distinct series (or two series 
in AD) as in general there are several different volumes Njj corresponding to a given 
Nq 1 ° ,x . One of those series is an ensemble of stacked spheres, a particular simple 
class of triangulations which can be enumerated exactly. This we do in Section 3. 

The other two series can be considered as "almost" stacked spheres, i.e. with 
one, or two, defects respectively. Heuristic arguments, supported by numerical ev- 
idence, show that although the leading asymptotic behavior of the three series is 
identical, the effect of inserting a defect increases the exponent 7 by +1 (or +2) 
compared to stacked spheres. This is discussed in Section 4. 
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In the second part of the paper, Section 5, we study the model Eq. (g) with a 
modified measure, 

W max (N D ,P) = — H qf. (6) 

Here qi is the order of a vertex i, i.e. the number of simplices incident to that vertex. 
As the measure is modified, varying /3, we observe a soft geometric phase transition 
from the elongated phase to a crinkled phase for (3 < —1, in agreement with the 
phase diagram Figure 1. Moreover, a scaling analysis of the transition is consistent 
with a third order transition (a = — 1). 

In the crinkled phase the fractal structure is dominated by a gas of sub-singular 
vertices, i.e. vertices whose order grows sub-linearly with the volume. The intrinsic 
fractal structure is characterized by a set of critical exponents: (i) A negative string 
susceptibility exponent 7, which decreases with j3. (ii) A Hausdorff dimension 
which is either dn = 00, or d' H ~ 2 (as j3 — > —00), depending on whether it is 
measured on the triangulation itself or its dual graph, (iii) A spectral dimension, 
measured on the dual graph, which increases smoothly from a branched polymer 
value d s = 4/3, at (3 = 0, to d s « 3/2 at /3 = -4. 

We have collected some of the more technical issues in two appendixes: In 
Appendix A we list the SCE of the three different series, and in Appendix B we de- 
scribe some complications in extracting 7 from measurements of minbu distributions 
on maximal triangulations. 



2 Maximal triangulations 



A maximal triangulation is a triangulation which for given volume Nd has the 
maximal allowed number of vertices Nq. For a combinatorial triangulation this 
implies: 
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(7) 



Here [^J denotes the floor function, i.e. the biggest integer not greater than x. From 
this definition it follows that in general more then one volume will correspond to a 
given maximal vertex number. In 3D this defines three different series, labeled by 
the (No, N3) -pairs: 



S° (N Q ,3N -W), S 1 (N ,3N Q -9), 



and 



S 2 (N ,3N -8). (8) 



The first series S° has the smallest volume for a given number of vertices and we 
will call it the minimal series. 
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Figure 2: The (p, q)-moves in three dimensions: (a) Inserting a vertex and its inverse, 
deleting a vertex. (6) Replacing a triangle by an link and vice verse. 



Likewise, in AD there are in principle four different series. However, as only 
combinatorial four-triangulations of even volume are allowed this reduces to only 
two: 



S° (jVo,4jV -18) and S 1 (N ,4N - 17). 



(9) 



In the space of all triangulations, the different series are related through a 
set of (topology preserving) geometric changes, the (p, g) -moves. In a (p, q)-move, 
where p = D + 1 — q, a (q — l)-sub-simplex is replaced by its "dual" (p — l)-sub- 
simplex, provided no manifold constraint is violated. These moves are used in the 
MC simulations and are ergodic for D < 4, i.e. any two triangulations are related 
through a finite sequence of moves || . 

In three dimensions there are two sets of moves: Move f 1, 41 inserts a vertex 
into a tetrahedra and its inverse, move (4, l) , delet es a vertex of order four. In move 
(2,3) a triangle separating two adjacent tetrahedra is replaced by a link connecting 
the opposite vertices. This is shown in Figure g These moves change the number 
of different D-simplices in the following way: 



/N \ 

N 2 

w 



(M) 

(4,1) 



fN + 1\ 

iVi +4 
iV 2 + 6 
V^s + 3/ 



and 



f N °\ 
W 



(3,2) 



JVi + 1 
N 2 + 2 
\N 3 + I J 



(10) 



Note that only the first set of moves changes both the volume and the number of 
vertices, and when applied to a maximal triangulation we stay within the particular 
series, i.e. the different series are internally connected via this move. 

The different series are connected, for fixed No, by the (2,3) and the (3,2)- 
moves. However, as it is not possible to apply move (3, 2) to the series S° (as 
triangulations in this series have the minimal volume for a given No), the only 
connection which the maximal ensemble has with the rest of the triangulation space 
is via move (2, 3) applied to the series S 2 . The relations between the different series, 
via the (p, g)-moves, are shown in Figure 3. 
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Figure 3: The relations between the different series of the 3-D maximal ensemble, S°, S 1 
and S 2 , via the (p, g)-moves. 

2.1 Stacked spheres 

A D-dimensional stacked sphere is a triangulation constructed by successively 
gluing together smallest D-spheres made out of D + 1, Z)-dimensional regular 



simplexes[10]. Gluing is defined as cutting away one .D-simplex in each triangulation 
and joining the two triangulations by the resulting boundary. Gluing a D-sphere 
corresponds to applying the (1, D + l)-move. 

From the definitions of stacked spheres and the maximal ensemble it follows 
that stacked spheres belong to the minimal series S°. The inverse statement, al- 
though less trivial, is also true. This was proven, in three and four dimensions, in 
two theorems due to Walkup |Q : 

Theorem 1 For any combinatorial triangulation of a 3-sphere the inequality 

N 3 > 3N - 10 (11) 
holds with equality if and only if it is a stacked sphere. 



Theorem 2 For any combinatorial triangulation of a ^sphere the inequality 

N A > 4N - 18 (12) 
holds with equality if and only if it is a stacked sphere. 

From those theorems follows a non-trivial statement: The move (l,D + l), 
and it inverse (D + l, l) , is ergodic when applied to the ensemble of stacked spheres. 
In contrast, for the non- minimal series this statement is not true. There exist 
triangulations, T G S k , k > 0, that cannot be reduced to the minimal configuration 
of the corresponding series by a repeated application of the move (-D +1, l). This 
we have verified numerically. 

2 The ensemble of stacked spheres is as well defined in two dimensions as for D > 2. However, in 
2D this does not define a maximal ensemble as ./V2 = 2Nq — 4 for every (spherical) triangulation. 
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3 Enumeration of stacked spheres 



As stacked spheres have a simple tree-like structure, it is possible to calculate the 
partition function restricted to this ensemble explicitly. In this section we present 
algorithms for doing so, both for labeled triangulations (which enter the MC simula- 
tions) and for the number of distinct or unlabeled triangulations. We also show how 
to generalize the algorithm for counting labeled triangulations to include a modified 
measure Eq. (||). 

All the results in this section are based on Ref. [12] which deals with the enu- 



meration of simplicial clusters. A (D + 1) -dimensional simplicial cluster is a simpli- 
cial complex obtained by gluing together (D+l)-dimensional regular simplices along 
their D-dimensional faces. Comparing this to the definition of stacked spheres one 
sees that a .D-dimensional stacked sphere is the boundary of a (D + l)-dimensional 
simplicial cluster where the outer faces of the simplicial cluster correspond to the 
simplices of the stacked sphere. So number of (D+l)-dimensional clusters with n 
(D + l)-simplices corresponds to the number of stacked spheres with Nq = n + D + 1 
vertices, i.e. Nd = Dn + 2 simplices. 

3.1 Labeled triangulations 

The basic "building block" is the expression for the number of (D + l)-dimensional 
simplicial clusters build out of n (D + l)-simplices rooted at a marked outer face. 
Denoting this number by en+i,n we have the following recursive relation [O], 

(13) 

niH \-n D+1 =n—l 

Introducing the counting series, 

oo 

E D+ i(t) = e D+hnt n = l + t(E D+1 (t)) D+ \ (14) 

n=0 

from Lagrange's inversion formula (see for example |]l3| page 147) it follows that 

1 f(D + l)n\ 

e D +l,n = n , : • (15) 

nD + 1 V n J 

To calculate the number of labeled simplicial clusters we note that once the 
outer face has been labeled, the remaining Nq—D vertices can be labeled in (Nq — D)\ 
distinct ways. The four labels on the marked face can be chosen in one of (p) 
ways. Including the factor from Eq. (||), a normalization factor (D + 1)! 
(as used in Ref. M), and dividing by Nd to adjust for the fact that we count 
rooted configurations, we arrive at the micro-canonical partition function for stacked 
spheres, 

Ws(N D ) = 2+1 e D+1 ^ . (16) 
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The first 20 terms in the above series are listed in Tables 2 and 3 in Appendix A for 
D = 3 and 4. This corresponds to the enumeration of all stacked spheres of volume 
N 3 < 62 and iV 4 < 82, respectively. 

To investigate the leading behavior of the above series, we perform an asymp- 
totic expansion of Eq. ([u]) in powers of N^ 1 . For D = 3 and 4 this gives: 

jV 3 -2 

«r 10 /256\^~ / 83 1 11305 1 

Wsm = {*r) [ 1 + tsn; + msNi < 17 > 

2109275 1 ^ 503051857 1 30688689031 1 
+ 663552 iVf + 127401984 ivf + 6115295232 ivf + "' 

_ iV 4 -2 

T,r/»rN /3125\^/ 33 1 1729 1 

WsiNl> = ( » ) V 1 + «W t + WjvJ < 18 > 

201159 1 19060503 1 556062507 1 



' 80000 iVf + 6400000 iVf + 128000000 Af + "' 

In both cases, the leading asymptotic behavior is of the form Eq. (|j) with 7 = 1/2. 
Moreover, for both series the corrections to the leading order are analytic in N^ 1 , 
i.e. there are no confluent singularities. This behavior is essential for applying the 
particular variant of the ratio method we use to analyze the series in next section. 

In contrast, the enumeration of the unlabeled triangulations — the number of 
distinct graphs — is much harder as the symmetry factor of each triangulation has 
to be canceled. This is the main subject of Ref. [12] which provides an algorithm to 
construct the generating function 

00 

h D +l(t) = h D+l,nt n , (19) 

n=0 

where Hd+Iu ls the number of distinct (-D+1) -dimensional simplicial clusters build 
out of n simplices. We list the functions h^(t) and h${t) in Appendix A together 
with the first 20 terms (Tables 2 and 3). It is interesting to note that although the 
series of unlabeled triangulations has the same asymptotic behavior as the labeled 
one, i.e. the same fj, c and 7, the sub-leading corrections are not analytic in this 
case. This is evident from an analysis of the singularities of the functions h^(t) and 
h$(t). As a consequence, for unlabeled triangulations it is is impossible to extract 
the leading behavior using a series extrapolation method such as the ratio method 



(as is demonstrated in Figure 13). 



3.2 A modified measure 

It is possible to extend the above algorithm for enumerating labeled stacked spheres 
to accommodate the measure term q 13 . The additional complication is to keep track 
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of the order of each of the vertices on the marked face. Here we present a recursive 
algorithm to do this. However, for simplicity we consider only 3D stacked spheres, 
i.e. 4D simplicial clusters. 

We denote by e(q; n, /3) the sum over all triangulations rooted at a marked 
outer face, where q = qi, ■ ■ ■ ,04 is the order of the vertices at the root. Each 
triangulation is weighted by a factor 

No 

n of, (20) 

i=5 

where the product is over all vertices not belonging to the marked face. Thus the 
number of rooted simplicial clusters becomes 

e(q;n,/3) = ^ ^ e(a;ni,/3) e(b;n 2 ,(3) e(c;n 3 ,/3) e(d;n 4 ,/3) 

»i+«2+"3+™4 a,b,c,d „ 

=n - 1 x (a 4 + 6 4 + c 4 + fi 4 -4) p (21) 

x ^ai+6i+ci-2,iji ddi+a 2 +b 2 -2,q 2 ^c 2 +d 2 +a 3 -2,qa ^b 3 +c a +d 3 -2,q4, > 

with the normalization 



e(2,2,2,2;0,/3) = 1. (22) 
We have to sum over all possible vertex orders of the marked face, 

e(n,f3) = (^i 92 Q3 Qa) P e(q; n, 0) , (23) 

q 

to obtain the micro-canonical partition function 

W S (N 3 ,(3) = 1 JL e (*p,0), (24) 

where = 4 5 ^ is the weight of the smallest triangulation. 

Although a solution of Eq. ([H|) is not known in a closed form, it is possible 
to calculate the series Ws(N 3 , (3) recursively. A simple program written in Fortran 
90 using quadruple precision takes about 7 hours of CPU time on a 600MHz Alpha 
station to compute Wg(iV*3,/3) for ./V3 < 53 (n < 17). Each additional term in the 
series requires four times the amount of computation as the previous one. 

Compared to a method of directly identifying distinct triangulations in MC 
simulations, discussed in next section, using the algorithm presented above it is 
possible to calculate more terms in the SCE. However, method presented above is 
limited to the enumeration of stacked spheres, i.e. to a sub-class of the maximal 
ensemble, and moreover the series has to be calculated separately for each value of 
P. 
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4 Pure maximal triangulations 



We will start our exploration of the maximal ensemble with the pure (unmodified) 
case P = 0. As stated, the partition function Eq. (||) effectively separates into three 
distinct series, S°, S 1 and S 2 , corresponding to regions of the triangulation-space 
only loosely connected via the (p, g)-moves. In this section we would like to explore 
in more details the difference between those series, in particular their asymptotic 
behavior. 

The following heuristic argument gives some idea of what asymptotic behavior 
we can expect. The series S 1 and S 2 , which are close to but do not saturate the upper 
bound in Walkup's theorem Eq (|TT|), can be considered as stacked spheres with one, 
or two, defects, respectively. A defect is defined as an one application of the move 
(2,3), i.e. replacing a triangle by a link. For a sufficiently large triangulation one 
naively expects that the number of ways such a defect can be inserted should grow 
linearly with the volume (inserting defects into two distinct non-symmetric stacked 
spheres does not lead to identical triangulations). This implies that although the 
leading exponential behavior of the three series might be the same, the sub-leading 
behavior is should be multiplied by A3 (or A|), hence 7(5' 1 ) = 1.5 (or j(S 2 ) = 2.5), 
with 7 defined in Eq. (Q). 

To investigate the asymptotic behavior, and to verify the statement above, 
we have performed several numerical experiments. In Figure |] we show the micro- 
canonical partition function for the three series measured numerically for A3 < 300. 
We have divided out the leading behavior of the partition functions, i.e. we plot 

S k = c k e-M*0^3 n 3 3 -^ s k , (25) 

where Ck is some convenient normalization and we assume the exact values 

H c (k) = - log f —J and 7 (fc) = ~ + k. (26) 

The behavior in Figure |i| is indeed compatible with this assumption. If the series 
had different /J, c (k) the curves should approach a straight line with a non-zero slope. 
Likewise, a different value of j(k) would manifest itself in a very slow convergence. 
Neither behavior is observed, all three curves in Figure || converge rapidly on a 
horizontal line. It should be noted, however, that the two non-minimal series appear 
to have larger finite-size corrections. 

4.1 The strong- coupling expansion 

A more precise determination of the critical parameters, p c {k) and ^(k), comes 
from the SCE of the micro-canonical partition function Eq. (||). In Section 3 we 
presented a algorithm for calculating the SCE restricted to the minimal series S . 
A corresponding recursive algorithm for constructing the two non-minimal series is, 
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Figure 4: The partition functions S of the three distinct series of the maximal ensemble 
in D = 3. The leading behavior, e Mo ^ fe ^ 3 N^ k ^~ 3 , has been divided out and the curves 
are normalized with an estimate of the value at iV3 = oo. The symbols just identify the 
different series. 

however, not available. In fact, it is not clear that such algorithm exist. Unlike 
5°, the non-minimal series cannot be constructed from some (different) minimal 
configuration by a successive application of move (1,4). 

We thus resort to an alternative method to enumerate the two non-minimal 
series. This method, based on a direct identification of all distinct triangulations 
in MC simulations, has previously been applied in four dimensions to identify all 
distinct triangulations of volume N± < 38 |B|. The two main ingrediences of the 
method are: (i) A sufficiently complicated hash function f(T) that uniquely identi- 
fies combinatorially distinct triangulations. (ii) A calculation of the corresponding 
symmetry factors, Ct, by an explicit permutation of the vertex labels. 

The SCE's for Si and 52 are shown in Table 4 in Appendix A. As the two 
non- minimal series grow much faster than the minimal one, this method is limited 
to volume N% < 36. This corresponds to 10 and 8 terms respectively. To calculate 
the next term would require the identification, and storing, of ~ 3 x 10 7 distinct 
triangulations — a formidable task with the computer resources available. 

To analyze the SCE and to extract the asymptotic behavior we apply an ap- 
propriate series extrapolation method. For these particular series the most powerful 
method is a variant of the ratio method, a Neville- Aitken extrapolation, in which 
successive correction terms to the partition function are eliminated by a suitable 
sequence extrapolation technique [Q. For the minimal series this method yields 
impressively accurate estimates of the critical parameters. Applied to the first 20 
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terms in Table 2 we get 



H C (0) = 0.749780192825077800383(15) , 
7(0) = 0.500000000000000000002(34) , 

compared to the exact values 7 = | and 

j» c (0) = 0.74978019282507780038404 . . . 

The quoted errors indicate how the estimate change as the last term in the series is 
excluded. 

For the non-minimal series the finite-size corrections are larger and, in ad- 
dition, we have less terms available. The estimates of the critical parameters are 
correspondingly less accurate, especially for S 2 : 

/i c (l) = 0.749765(17) 1 x 
7(1) = 1.497(18) J ' 

fi e (2) = 0.74940(49) \ 2 
7(2) = 2.24(30) J 

These values are though consistent with the conjecture that fi c is the same for all 
the three series, but that j(k) = | + k. 

The same behavior is observed for the two series of the maximal ensemble 
in AD. The first 20 terms in Table 3 yield equally impressive agreement with the 
exact values: \i c = log ^3125/256 ~ 0.625503... and 7 = i. And the first 7 
terms in the non-minimal series S , given in Ref. ||], give /x c (l) = 0.6245(17) and 
7(1) = 1.64(35). 

Some final comments regarding the non-minimal series. Although an exact 
enumeration algorithm is not available it is possible to obtain a lower bound on the 
series S* 1 and S 2 by assuming that all non-minimal triangulations can be constructed 
by move (l,4) from the appropriate minimal configuration. The resulting SCE has, 
for both series^, the same [i c and 7 = | . 

4.2 Distributions of baby universes 

The different asymptotic behavior for the three series is very pronounced when one 
considers the distribution of minimal neck baby universes (minbu's) measured on 
maximal triangulations. A minbu is a part of a triangulation connected to the rest 
via a minimal neck — on a three dimensional combinatorial triangulation a minimal 
neck consists of four triangles glued along their edges to form a closed surface (a 
tetrahedra) . By counting in how many ways a "baby" of size b can be connected to 

3 Note that for the series S 2 this assumes that the two defects are connected, thus excludes 
triangulations with two randomly distributed defects. This results in a lower estimate 7=5, 
instead of 7 = | which we expect for S 2 . 
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Figure 5: The minbu distributions, B 0,a '° and S 1 ' 1 ' 2 , measured on triangulations of size 
N Q = 170 (i.e. N 3 = 500 or 502). In the figure we plot y = \ogB vs. x = log (6 (1 - b/N 3 )), 
and the dashed lines are fits to Eq. (Ejj). 



a "mother" of size N$ — b + 2, the distribution of minbu's can be expressed in terms 



of the canonical partition functions of the corresponding volumes [15|. 

For a maximal triangulation the situation is more complicated as the baby, 
the mother, and the "universe" can belong to any one of the different series S k . It 
is not even obvious that if a maximal triangulation is split along a minimal neck the 
two parts belong to the maximal ensemble at all. This is though easy to verify^]. 
In general the distribution of baby universes measured on maximal triangulations 
becomes 

B k B ,k M ,k u(b) _ (b) bS k Hb) (N 3 -b + 2) S ""(N 3 -b + 2) 

~ 6 7(fcfl) ~ 2 (JV 3 - b + 2) 7(fcM) ~ 2 . (27) 

The indices fee, 1cm, and kjj refer to the baby, the mother, and the universe, the 
factors b and N3 — b + 2 count the possible locations of the neck, and Cn 3 (b) is a 
contact term due to interactions on the neck. The contact term is discussed further 
in Appendix B; for pure simplicial gravity this is just a trivial constant. In the last 
step we have used the expected asymptotic behavior of the series Eq. (|j). 

In principle, Eq. ( p7|) defines nine different distributions, however the addi- 
tional restriction, ks + fcju < 2, reduces the number of possible combinations to 

4 The converse is however not true. If we glue together two maximal triangulations, Ti £ S kl 
and T2 G S* 2 , the result only belongs to the maximal ensemble provided k\ + ki < 2. 
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six. Two of those, £>°> ' and B ,1 ' 2 , both corresponding to ks = ku-, are plotted 
in Figure |f| Both distributions are measured on triangulations with iVo = 170, but 
on different volumes, N3 = 500, and 502. The former distribution, B 0,0 ' , which is 
composed solely of the minimal series, has remarkably small finite-size effects. A 
fit to Eq. (p7|), including minbu's b > 10, gives 7(0) = 0.495(4). In contrast, for an 
acceptable fit to the distribution B > 1 ' 2 a lower cut-off, b > 50, had to be imposed. 
The result, 7(1) = 1.44(4), is consistent with our previous estimate. 

The remaining four distributions in Eq. (27) correspond to a baby and a 
mother belonging to different series. However, the mixture of different 7's and 
large finite-size effects makes a reliable estimate of "f(k) very difficult from these 
distributions. 



5 A modified measure 

We now turn to the maximal ensemble modified by varying the measure Eq. (^) for 
—4 < (3 < 2. For most part we restricted our investigation to the ensemble of stacked 
spheres where we expect less finite-size effects. This provides the additional benefit 
that the MC simulations are much simplified as only move (l, 4), and its inverse, 
are required for the updating process to be ergodic. However, we also did extensive 
simulations of the full model (including finite values of kq) to verify that this did 
not introduce any bias in the measurement process. Typically we collected about 
10 to 50 thousand independent measurements at each value of (5 and on different 
volumes N 3 = 125 to 8000. 

5.1 Evidence for a phase transition 

As the fluctuations in the vertex orders are enhanced by modifying the measure, 
i.e. by decreasing /3, the internal geometry changes. The polymerization of the 
triangulations is suppressed and the model enters a new phase akin to the crinkled 
phase observed in four dimensions ||. The nature of the internal geometry in the 
crinkled phase is discussed in Section 5.2, here we present what evidence we have 
for the existence of a continuous phase transition at a critical value (3 C . 

In Figure |6]a we show the fluctuations in the energy density, Cy = ^((e 2 ) — 
(e) 2 ) , as we vary [5. This is for A3 = 125 to 8000 and we define e = N^ 1 Y^i log(gi). 
We observe on all volumes a peak in Cy consistent with the existence of a phase 
transitions at (3 ~ —1. The peak value appears to saturate as the volume is increased; 
this makes it difficult to distinguish between a soft continuous phase transition and a 
cross-over behavior. The scaling of the peak value is, however, consistent with a third 
order phase transition. This is shown in Figure ^ where we fit the maximal value to 
the scaling behavior: C™ ax ps a + bN^ ,udH . The fit yields a/vd H = -0.34(4), with 
a xV(d.o.f.) « 0.8. Assuming hyper-scaling, a = 2 — udn, is valid this implies a 
specific heat exponent a = —1.03(9). 
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Figure 6: Evidence for a phase transitions as the measure is modified varying (3. (Top) 
The fluctuations in the energy density Cy — -^((e 2 ) — (e) 2 ), where e = N$ z2i^°s(Qi)- 
(Bottom) The fluctuations in the maximal vertex order, \q = ^3\( e( lo) ~ ( e )( ( ?o))- 



Additional evidence for a transition is provided by other geometric observ- 
ables such as the maximal vertex order qo whose value jumps at the transition. As 
discussed in the Introduction, and is quantified in the Section 5.2, the crinkled phase 
is characterized by a gas of sub-singular vertices, i.e. vertices whose local volume 
grows like N% , x < 1. As a consequence both the fluctuations in the maximal vertex 
order, x qo = N 3 ((($) - (qo) 2 ), and its energy derivative, D ?0 = N 3 ((q e) - (qo)(e)), 
diverge at the transition. An example of the former is shown in Figure ^b. In fact, 
the susceptibility Xq appears to diverge in the whole crinkled phase in contrast to 
its behavior in the elongated phase. 

Note that some of the data in Figure |6| are measured on non-minimal triangu- 
lations (e.g. A3 = 1000 and 4000). This shows that the phase transition is a generic 
behavior of the maximal ensemble, not just of the minimal series. The same is true 
for the nature of the crinkled phase explored in the next section. 
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Figure 7: The volume scaling of the peak in the energy fluctuations, C™ ax . The curve is 
a fit to the scaling form: a + bN£; the optimal fit yields x — a/vdn = —0.34(4) with a 
x7(d.o.f.) a 0.8 

5.2 The crinkled phase 

To investigate the nature of the internal geometry in the crinkled phase, and to 
compare it to the elongated phase, we have measured several aspects of the fractal 
structure: the distribution of vertex orders ir(q), the string susceptibility exponent 
7, and both the fractal and spectral dimensions dn and d s . 

5.2.1 The vertex order distribution 

The shape of the vertex order distribution Tr(q) is a good indicator for the nature 
of the internal geometry in the different phases. In the elongated phase the vertex 
orders are more or less equdistributed and the maximal vertex order grows logarith- 
mically with the system size, go ~ log -/V3. In the crinkled phase, on the other hand, 
vertices of large order are favored and the maximal vertex order grows sub-linearly, 
qo ~ N^. Estimates of a, from finite-size scaling including volume N$ = 125 to 
8000, indicate that a ~ 0.8 for sufficiently large negative (5. 

Examples of the vertex order distribution in the two phases are shown in 
Figure ||for j3 = and = —4. In the crinkled phase the distribution 7r(g) developes 
a "bump" in the tail, corresponding to a condensation of vertices of large order. 
The bump persists as the volume is increased, i.e. it captures a finite fraction of the 
volume as -/V3 — > 00. This suggests that the crinkled phase can be characterized as 
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Figure 8: A log-log plot of the (normalized) vertex order distribution, n(q), both in the 
elongated phase, (3 = (Top), and in the crinkled phase, /3 = — 4 (Bottom), and measured 
on different volumes. 



a gas of sub-singular vertices.^ 

5.2.2 The string susceptibility exponent 7 

To measure the string susceptibility exponent for various values of /3, and restricted 
to the minimal series, we have used several methods: 

(a) The pseudo-critical cosmological constant, fJ, c (Ns) 

Assuming the asymptotic behavior Eq. (||), a saddle-point approximation of the 
grand-canonical partition function gives: 

Hc(N 3 ) « fi c + . (28) 



Similar behavior, although more pronounced, is observed in the crumpled phase in AD when the 
ensemble of triangulations in Eq. (^) includes degenerate triangulations However, in 4D the 

bump is disconnected from the rest of the distribution — this justifies the nomenclature singular 
structure. 
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Table 1: Estimates of the exponent 7, for /? = — 1, —2, —3, and —4; (a) from a fit to /z c (7V 3 ), 
Eq. (p8|), (with [x 2 /(d.o.f.)]); (b) from the minbu distribution, Eq. (p7|), both including and 
excluding the contact term, B 0,0 ' and i5°'°'°; and (c) from the first 17 terms in the SCE. 





(J=-l 


P = -2 


p= -3 


= -4 


(a) 




-0.37(4) [1.5] 


-2.56(3) [0.8] 


-5.25(2) [0.6] 


-8.11(4) [0.7] 


(b) 


£0,0,0(6) 


0.25(7) 


-0.20(3) 


-0.37(4) 


-0.52(8) 




B°A°(&) 


-0.14(5) 


-2.58(7) 


-6.5(2) 


-10.5(8) 


(c) 


SCE 


-0.076(23) 


-2.401(17) 


-5.484(42) 


-8.455(23) 



By measuring fJ, c (N^) with high accuracy on volumes N$ ~ 100 to 2400, and for 



(3 = — 1, -2, -3, and -4, we determined 7 by a fit to Eq. (28). 

(b) The minbu distribution, B°'°'°(b) 

We measured the minbu distribution on volume N3 = 2000 and for the same 
values of (3. However, in contrast to the elongated phase, in the crinkled phase 
the contact term Cjv 3 (&) has a non-trivial dependence on b and has to be elim- 
inated prior to extracting 7 from the (corrected) distribution £>°>°'° by a fit to 
Eq. (f27|). 

(c) The strong-coupling expansion I 

Using the enumeration methods for a non-zero /3, presented in Section 3.2, we 
calculated the first 17 terms of the SCE of stacked spheres for few values of [3. 
From this SCE we extracted 7 using the ratio method. 

(d) The strong-coupling expansion II 

By explicitly identifying all distinct triangulations of volume N3 < 38, and 
calculating the corresponding measure term, we calculate the first 12 terms in 
the SCE for an arbitrary value of (3. 

The estimates of 7 are listed in Table 1 and plotted in Figure |9[ Qualitatively, the 
agreement between the different methods is very good. Only the estimate from the 
(corrected) minbu distribution B 0,0,0 deviates for large negative [3. This is under- 
standable as for large negative (3 the tail of the minbu distribution is suppressed 
and is difficult to measure accurately. 

Note that if we do not correct the minbu distribution for the contact term 
Cat 3 (6), for (3 < a fit to £? '°> yields a very different, and unreliable, estimate of 
7. This estimate in included in Table 1 and Figure ^. Why the contact term is 
so important for maximal triangulations in the crinkled phase, and how it can be 
eliminated, is discussed in Appendix B. 
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The agreement between the different methods strongly supports the assump- 
tion that the asymptotic behavior of the canonical partition function, Eq. (||), is as 
valid in the crinkled phase as in the elongated phase. There are, however, indica- 
tions from the SCE that this might not be true at the phase transition. This we 
show in Figure [h]. In the upper part of the figure we show the variation of 7 with 
/3, estimated from the first n terms in the SCE, for n = 6, 8, 10 and 12. In the lower 
part we plot an estimate of the convergence of the ratio method, 

A 7 n = - i-L 1 L . (29) 

' 2 |7 n | + | 7 n_1 | 

In the elongated phase and in the crinkled phase the estimates of 7 converge. How- 
ever, both for « — 1, and for 3> 0, the method fails to converge. It is possible 
that at the phase transition there are some different non-analytic finite- volume cor- 
rections to the asymptotic behavior Eq. (Q). There are known examples of such 
behavior: For one scalar field (c = 1) coupled to 2D-gravity there are logarithmic 
corrections to the free energy fl7|| . 

Why the method appears to fail for (3 3> is less clear. As there is no phase 
transition in this region, we would expect Eq. ([|) to be equally valid for (3 > as for 
(3 = 0. We note, however, that in general the magnitude of the finite-size corrections 
increases drastically for large positive f3. 
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All the above estimates of 7 are for the minimal series. We expect, however, 
that the arguments presented in Section 4 are also valid for non-zero values of /?, 
and that for a given /3, 7 is correspondingly larger for S 1 and S 2 than for S°. The 
determination of 7 is though more difficult in this case: we have no algorithm for 
enumerating the non-minimal series; the minbu distributions are more complicated 
(as is discussed in Section 4.2); and as the MC simulations are not ergodic restricted 
to either of the non-minimal series, we cannot determine ^i c (N^) numerically. 

We are thus left only with method (d). As for S°, we have identified all distinct 
non- minimal triangulations of volume N3 < 36; re- weighted with the measure this 
yields the variations of j(k) with 0. For example, for f3 = —2 we get 7(5' 1 ) = 
-1.2508(36) and 7(5 2 ) = -0.15(33), compared to 7(5°) = -2.401(17). Similar 
results are obtained for other values of (3. In all cases the estimates of jik) are 
consistent with our expectations of how the three series behave. 
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5.2.3 The Hausdorff and the spectral dimension 



We have estimated both the Hausdorff and the spectral dimension of the maximal 
ensemble in the elongated and in the crinkled phase. On smooth regular manifolds 
those two definitions of dimensionality coincided, however on highly fractal random 
geometries, likes the ones that dominate the partition function Eq. (H), they are in 
general different. 

The Hausdorff dimension, dn, is defined from volume of space within a ball 
of geodesic radius r from a marked point: V(r) ~ r dH . We have measured this 
both on the dual graph, from the simplex-simplex distribution sn 3 (t), and on the 
triangulation itself, from the vertex- vertex distribution vj^ (r). These distributions 
count the number of simplices (vertices) within a geodesic distance r from a marked 
simplex (vertex). Assuming that the only relevant length-scale in the model is 
defined by x = N s l ^ dn , general scaling arguments imply 

SN 3 (r) = N l f l/dH F{x) , (30) 

and likewise for VN (r). By comparing distributions measured on different volumes, 
Eq. (|30|) provides an estimate of djj. 

In principle, the Hausdorff dimension defined on the dual lattice, d' H , need not 
to agree with the one defined by Vjv (r), du- However, they do so in most models 
of dynamical triangulations which have a continuum interpretation as theories of 
extended manifolds. This is also the case in the elongated phase where we expect 
du = 2, characteristic of generic branched polymers. We measured the above distri- 
butions on triangulations of volume N3 = 125, 500, 2000 and 8000, and determined 



dn from the scaling Eq. (p0|). The estimates are plotted in Figure 11. 

But, not unexpected, as is decreased and we pass through the phase tran- 
sition the Hausdorff dimension changes. What is surprising though is that in the 
crinkled phase the estimates of d' H and du differ substantially. The former, mea- 
sured on the dual graph, jumps to d' H 3 at the transition, then decreases with @, 
possible to d' H ~ 2 as (3 — > —00. In contrast, the vertex-vertex distribution yields a 
very large estimate, dn ~ 10 for (3 < —2. Such large estimates on a finite volume 
suggest that the actual value is infinite. 

This discrepancy is much too large to be explained by an uncertainty in the 
determination. It is correct that in the crinkled phase the vertex- vertex distribution 
has an extends of only few geodesic steps[] (all vertexes are very "close"), even on 
volume 8000. This makes the determination of du correspondingly difficult. How- 
ever, this is just a signal of an infinite fractal dimensions.. The determination of 
d' H from s(r) is, on the other hand, very reliable and, if anything, the estimates ap- 
pear to decrease with volume. We show examples of both distributions in Figure 12 



6 As a consequence, only the two largest volumes, corresponding to iVo = 670 and 2670, are used 
in determining du- And the volume scaling was not sufficiently good for a reliable error estimate. 



21 



5 - 



4 - 



3 - 



2 - 



d 



s 

1.5 
1.4 
1.3 
1.2 



O {125,500} 

□ {500,2000} 

O {2000,8000} 

• {670,2670} 



O 



O 



~* — 



-r - 

-3 



T — 1 — I — ' — I ■ I — ■ I ■ I 



d 



H 



-10 

- 8 

- 6 

- 4 

- 2 

- 







1 fi 








500 


® CD 


□ 


2000 


M I | g 






o 

i 


8000 








w 

CD 

O 



Figure 11: (Top) The Hausdorff dimension denned either by the (a) the simplex-simplex 
distribution s(r), or (b) the vertex- vertex distribution v(r), d' H (open symbols) or dn (filled 
symbols). Note that the scale for dn is double compared to d' H . Triangulations of size 
{N3, 47V 3 } (or {No, 4iV }) are used in the scaling Eq. (j3p|). (Bottom) The spectral dimension 
d s measured on the dual graph. In both plots the horizontal lines is the generic branched 
polymer value. 
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Figure 12: (Left) The simplex-simplex and (Right) the vertex-vertex distributions, s(r) and 
v(r), for (3 — —4 and N3 = 500, 2000 and 8000. Both distributions are scaled according to 
Eq. ( p0| ) and are plotted vs. the scaling variable x = r/N^ dfI . Notice, however, that different 
fractal dimensions, d' H = 2.36 and djj — 12, are used in scaling the two distributions. The 
solid line is a polynomial interpolations. 



for (3 = —4. Both distributions have been scaled according to Eq. (30) using the 
corresponding optimal values of dn- 

Finally, we also determined the spectral dimension d s from the return proba- 
bility of a random walker on the dual graph: p(t) ~ t~ ds / 2 . The time t is measured 
in units of jumps between neighboring simplexes, with a hopping probability 1/4. 
These estimates of d s , for volumes N3 = 500, 2000, and 8000, are included in Fig- 
ure 11. In the elongated phase we get d s ~ 4/3, the value for a generic branched 
polymer, whereas in the crinkled phase the spectral dimension increases and reaches 
d s ~ 3/2 for (3 = —4. Similar behavior is observed for the spectral dimension in 
AD ||, although in that case the estimate in the crinkled phase was somewhat 
larger (d s ~ 1.8 for (3 = —3.5). It is notable, however, that both in three and four 
dimensions the spectral dimension in the crinkled phase is less than two. 



6 Conclusions 

In this paper we have investigated an ensemble of maximal triangulations i.e. tri- 
angulations that, for a given volume are as close to the upper kinematic bound, 
Nq/Nd = as possible. This is motivated by the observation that this ensem- 

ble appears to reflect the properties of the weak-coupling limit of simplicial gravity, 
Kq — > 00. 

In the first part of the paper, we presented an algorithm which, in the un- 
modified case (3 = 0, yields an explicit enumeration of stacked spheres — a sub-class 
of the maximal ensemble. We also generalized this algorithm to include a mod- 
ified measure, which, for (3 7^ 0, allows a recursive construction of a SCE of the 
corresponding partition function Eq. @. 
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Using a combination of SCE and MC simulations, we investigated the prop- 
erties of the maximal ensemble. We found that it is composed out of three distinct 
series, all with the same leading asymptotic behavior, but with different leading 
corrections 7(fc). This can be understood as the maximal ensemble includes, in 
addition to stacked spheres, triangulations that are "almost" stacked spheres, i.e. 
stacked spheres with one or two defects. 

In the second part we demonstrated that with a suitable modified measure the 
maximal ensemble exhibits a transition to a crinkled phase, akin to what is observed 
in AD simplicial gravity. This agrees with the phase diagram Figure 1. The scaling 
properties of this transition are consistent with a third order transition, i.e. with a 
specific heat exponent a = — 1. In the crinkled phase the intrinsic fractal structure 
is charcterized by a string susceptibility exponent 7 < 0, and Hausdorff dimension 
dn = 00 or d' H 2, depending on whether it is measured on the triangulations, 
or its dual. This difference in the direct vs. dual fractal structure makes it unlikely 
that a thermodynamic limit exist in the crinkled phase, where a theory of extended 
manifolds emerges. The value of te spectral dimension measured on the dual graph 
agrees with the spectral dimension observed for the non-generic branched polymers. 
In the later case however exponent 7 is positive. In case of 7 < BP model predicts 
d s = 2 |1| 

We note the similarity between the above phase structure and what is observed 
in two dimensions if one varies the measure. There too is a soft continuous transition 
to a crumpled phase, for sufficiently negative (3, where the string susceptibility 
exponent is negative IS] and the internal geometry collapses. Indeed, the similarity 
between D = 2 and D > 2 extends further. In both cases, the same phase strucuture 
is observed when sufficent amount of matter is coupled to the model. This can be 
understood, qualitativly, as the discretized form of the matter fields includes to a 
leading order the measure Eq. @. 

The phase diagram Figure 1 can also be partially understood from an analysis 
of a mean-field model which is related to branched polymers [poj ] . Restricted to the 
maximal ensemble this model predicts a soft (third order or greater) phase transition 
between an elongated and a crinkled phase. In a language of branched polymers this 
corresponds to a transition from an elongated (or a generic) phase with 7 = 1/2 and 
du = 2, to a bush phase (non- generic) with 7 = 1-/3 and du = 00. The intrinsic 
crumpling in the bush phase is though much more pronounced than is observed in 
the crinkled phase of the maximal ensemble. 

To which extend do we expect the maximal ensemble to correctly reflect the 
properties of the whole weak-coupling phase, hq > Kq, of the model Eq. (Q). In 
principle, this may depend on how one takes the two limits: the thermodynamic 
limit N3 — > 00, and the weak-coupling limit kq — > 00. A priori, these two limits 
need not commute. However, in Ref. |l(| it is claimed that already for a finite 
kq > Kq, the partition function Eq. (0) is dominated by triangulations which are 
close to stacked spheres. 
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A Enumerations of the maximal ensemble 



In Tables 2 and 3 we list the first 20 terms in the SCE of both labeled and unlabeled 
staked spheres in three and four dimensions. The enumeration of unlabeled stacked 



spheres is given by the generating function Eq. (|19|), where: 



h^t) = ^(tE^t 3 ))- ^E 1 (tE 4 (t 3 )f+ ^E 2 (tE 4 (t 2 )) (31) 

- l E 2 (tEi{t 2 )) 2 + ^(t 2 E 4 (t 6 )) E 2 (t 3 £ 4 (t 6 )) + ±.E 4 (t) 

+ Ht 2 £ 4 (i 2 ) 4 + 1 -tE 1 (tE^)) 2 E i ^)+ \tE 4 (t±) 
+ h 2 E 4 (t*) 2 + \tE^) + \t 2 E^t 2 E^)) £ 4 (i 6 ) , 



h 5 (t) = -H+ l El (tE 5 (t 2 ) 2 )- ± El (tE 5 (t 2 ) 2 ) 2 + -^{tE^)) (32) 
-±E l (tE 5 (t i )) 2 + ±E 2 (tE 5 (t 3 ))- ±E 2 (tE 5 (t 3 )) 2 + ±E 3 (tE 5 (t 2 )) 

±E 3 (tE 5 (t 2 )) 2 + 1 tE 2 (t 2 E 5 (t 6 )) E 3 (t 3 E 5 (t 6 )) + ^E 5 (t) 
+ ±-t 2 E 2 (t 2 E 5 (t e )f E 3 (t 3 E,(t e )) 2 - ^E 5 (t) 2 + ^tE 5 (t) e 
+ i tE 3 (tE 5 (t 2 )) 4 E 5 (t 2 ) + 1 tE 1 (tE 5 (t 2 ) 2 ) 2 E 5 (t 2 f 
+ ^ 2 Mt 2 ) 5 + ^tE 2 (tE 5 (t 3 )) 3 E 5 (t 3 )+ ^tE 5 (t 3 ) 2 
+ ±tE 1 (tE 5 (t i )) 2 E 5 (t i )+ ±tE 1 (t 2 E 5 (t i f)E 5 (t i )+ ^tE 5 (t 2 ) 3 
+ lt 2 E 1 (t 2 E 5 (t i ) 2 )E 5 (t*)+ 1 -tE 5 (t 5 )+ ±t 2 E 5 (t 5 ) 2 
+ ltE 5 ^)+ U 2 E 2 (t 2 E 5 {t*)) 2 E,(f) + ±t 2 E 5 (t w ). 



The values in Tables 2 and 3 can be compared to Table 8 of Ref. [12]. Although for 
most volumes they agree there are few discrepancies. We have, however, an inde- 
pendent confirmation on our calculations; for Np < 38 we have explicitly identified 
all distinct triangulations. 

As discussed in Section 4.1, although the leading asymptotic behavior is iden- 
tical for the ensemble of labeled vs. unlabeled stacked spheres, for the latter the 
sub-leading corrections are non-analytic. This effect, which is demonstrated in 
Figure (|l3|), makes it virtually impossible to extract the asymptotic behavior of 
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Table 2: The first 20 terms in the SCE of the both labeled and unlabeled three dimensional 
stacked spheres, S°(N 3 ) and h^ n , respectively. This is for (3 — 0. 



n 


(N ,N 3 ) 




S°(N 3 ) 


1 


(5,5) 


1 


1 


2 


(6,8) 


1 


5/2 


3 


(7,H) 


1 


10 


4 


(8,14) 


3 


50 


5 


(9,17) 


7 


285 


6 


(10,20) 


30 


1771 


7 


(11,23) 


131 


11700 


8 


(12,26) 


795 


80910 


9 


(13 29) 


5152 


579700 


10 


(14,32) 


36800 


8544965/2 


11 


(15,35) 


272093 


32224114 


12 


(16,38) 


2077909 


247754390 


13 


(17,41) 


16176607 


1936016950 


14 


(18,44) 


127996683 


15340200750 


15 


(19,47) 


1025727646 


123020557800 


16 


(20,50) 


8310377720 


996993749142 


17 


(21,53) 


67967600763 


8155209690540 


18 


(22,56) 


560527576100 


67259885983090 


19 


(23,59) 


4656993996246 


558826638901000 


20 


(24,62) 


38949328897318 


4673871155753800 



unlabeled triangulations using series extrapolation methods. 

In Table 4 we list the first 10 and 8 terms of two non-minimal series, S 1 and 
S 2 , calculated by directly identifying all distinct triangulations in MC simulations. 

B The contact term 

As discussed in Section 5.2, for minbu's on maximal triangulations, the interactions 
at neck — the contact term — play a (unusually) big and non-trivial role in the 
crinkled phase. It is essential to take these interactions into account when estimating 
7 from a fit to the measured minbu distributions. In this appendix we analyze the 
behavior of the contact term and explain why it plays such an important role for 
this particular ensemble. 

In the nomenclature of Section 3 we can write the minbu distribution, Eq. (f27|), 
for (3 = and measured on a stacked sphere of size n, as 

B n (b) oc — '■ oc (3(n-&) + 2) (36 + 2) . , 
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Table 3: Same as Table 2, except for four dimensional stacked spheres. 



n 


(N ,N 4 ) 


^5,n 


S\N,) 


1 


(6,6) 


1 


l 


2 


(7,10) 


1 


3 


3 


(8,14) 


1 


15 


4 


(9,18) 


3 


95 


5 


(10,22) 


7 


690 


6 


(11,26) 


30 


5481 


7 


(12,30) 


142 


46376 


8 


(13,34) 


922 


411255 


9 


(14,38) 


6848 


3781635 


10 


(15,42) 


57994 


35791910 


11 


(16,46) 


525048 


346821930 


12 


(17,50) 


4999697 


3427001253 


13 


(18,54) 


49159506 


34425730640 


14 


(19,58) 


494873071 


350732771160 


15 


(20,62) 


5068890252 


3617153918640 


16 


(21,66) 


52632367550 


37703805776935 


17 


(22,70) 


552579655767 


396716804816265 


18 


(23,74) 


5855580019967 


4209161209968825 


19 


(24,78) 


62548009026369 


44993046668984145 


20 


(25,82) 


672818970206219 


484176486362971710 



However, for /3 ^ the situation is different. We must take into account the different 
possible distributions of vertex orders, p and q, on the marked simplex (the neck), 
both on the baby and on the mother. The distribution is then 

1 4 

B n (b; (3) oc — -— Y, J! & + « " 2 ^ ® e ^ b ' ® ■ ( 33 ) 



W 3 {n) 

y ' p,q i=i 



Defining 



e(q;n,/3) y-r p 

we can rewrite Eq. (p3|) as 



«(q;n,/3) = -^^LM > ( 34 ) 



oc (3(n-6) + 2)(36 + 2)- W) ^ (n - 6) 



W 3 (n) 

^ w(q.\n-b,(3)w(p;b,(3) f[ — + ^ g ^ 



/3 



B n (b;P)-C n (b;(3) , (35) 
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Figure 13: The successibe ratios, s^ 1 = S n /S n -i, for the sequence of labeled vs. unlabeled 
stacked spheres. The arrow indicates the correct asymptotic value z _1 = 256/27. For 
unlabeled triangulations the irregular behavior is due to non-analytic corrections. 



where we define C n [b; (3) as the contact term.. 

The following argument gives some idea of the relevance of the contact term 
in the different phases. If we assume that the distribution io(q) factorizes, 

4 

w(q;n,/3) = io(<fe; n, 0) , (36) 

then the distribution w(q; n) is related to the vertex order distribution n(q), i.e. the 
probability of finding a vertex with given order q, through 

/ m gir(q;n,f3) 
w{q;n,P) = — . (37) 

In the elongated phase the distribution w(q;n) decays exponentially, whereas in 
the crinkled phase vertices of large order are favored and the distribution has a 
power-law decay (Figure^). 

An exponential decay of w(q; n) effectively introduces a cut-off in the sum over 



vertex orders in Eq. (35). Hence, for sufficiently large triangulations the contact term 
becomes independent of the volume. In contrast, in the crinkled phase there is no 
such cut-off and we can expect a strong volume dependence. This is demonstrated by 



numerical measurements of the contact term, shown in Figure 14. In the elongated 
phase, [3 = 2, the contact term C n (2) is more or less constant for [3 > 50, whereas 
for [3 = —4 the contact term C n {— 4) depends strongly on b. 
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Table 4: The first few terms in the SCE of the 3D maximal ensemble for the two non- 



minimal series S 1 and S 2 . is the number of distinct (unlabeled) triangulations. 



N 


N 3 


Til 


S L (N 3 ) 


N 3 


"2 




6 


9 


1 


5/3 








7 


12 


2 


35/2 


13 


1 


15 


8 


15 


5 


152 


16 


8 


4205/16 


9 


18 


23 


3800/3 


19 


45 


3290 


10 


21 


124 


73370/7 


22 


385 


72535/2 


11 


24 


859 


348075/4 


25 


3435 


375750 


12 


27 


6518 


6544720/9 


28 


32710 


15043835/4 


13 


30 


52761 


6121632 


31 


312601 


36867985 


14 


33 


438954 


570959025/11 


34 


2995589 


713253255/2 


15 


36 


3717370 


886163135/2 









Although this effect can be expected for any ensemble modified by the mea- 
sure Eq. (|6]), the particular construction of maximal triangulations, by a repeated 
application of move (1,4), enhances the role of the contact term. For every time 
move (1,4) is applied one minimal neck is generated, hence for a maximal triangu- 
lation the number of minbu's is fixed. Thus the minbu distributions, measured on 
S°, is more constrained than measured on a general ensemble where triangulations 
with fewer minbu's are favored as (3 —* oo. 

How does all this affect the measurements of 7 ? The standard method for 
extracting 7 from a minbu distribution is to impose a lower cut-off bo on the minbu 
size included in the fit to Eq. (p7|). The cut-off is then increased until a stable 
estimate of 7(60) is obtained ]21]| . In the elongated phase, and in general for most 
models of simplicial gravity, this methods works fine. However, for stacked spheres 
this method fails in the crinkled phase. Unless the contact term is eliminated, i.e. 
divided out of the minbu distribution in the measurement process, a stable plateau 
in 7(60) is n °t observed. An example of this is shown in Figure [l4| for [3 = — 4. 
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Figure 14: (Left) The dependence of the contact term, C n (b; f3), on the minbu size b. This 
is for N3 = 2000 and both in the elongated phase, f3 = 2, and in the crinkled phase, (3 = —4: 
C(2) and C(— 4). (Right) The exponent 7, determined by a fit to the minbu distribution 
Eq. @ for 7V 3 = 2000 and (3 = -4. Results are shown varying a lower cut-off 60 in the fit, 
and both with and without the contact term divided out, 7 and 7 respectively. Note that 
in both plots the compared curves are plotted on very different scales. 
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